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Abstract 

This is a continuation of our study in [lOj concerning q-tori, i.e. low-dimensional tori in the phase 
space of nonlinear lattice models like the Fermi-Pasta- Ulam (FPU). In [10] we focused on the FPU-/3 
system, and we showed that the dynamical features of the q-tovi serve as an interpretational tool to 
understand phenomena of energy localization in the FPU space of linear normal modes (g-space). In 
the present paper we compute quasi-periodic g-tori solutions in the a model as well, using the method 
of Poincare - Lindstedt series. Then, we study the localization features of the q-tori solutions in both 
the a and (3 models. We consider various so-called 'initial excitations'. These are different choices of 
subsets of modes to which we give a non-zero oscillation amplitude at 0-th order of perturbation theory. 
We focus on two classes of initial excitations: i) ones describing low-frequency 'packets' of modes, and 
ii) excitations of a (small) set of modes with an arbitrary distribution in q-space. In the former case, 
we find formulae yielding an exponential profile of energy localization, following an analysis of the size 
of the leading order terms in the Poincare - Lindstedt series. In the latter case, we explain the observed 
localization patterns on the basis of a rigorous result concerning the propagation of non-zero terms in 
the Poincare - Lindstedt series from zero to subsequent orders. In both cases, the resulting localization 
patterns are found to agree to a large extend with the corresponding metastable states for so-called 
'FPU-trajectories', which have initial conditions close to q-torus solutions. However, we also discuss 
some differences between q-tori solutions and FPU-trajectories with similar localization profiles. Finally, 
we discuss the extensive (i.e. independent of the number of degrees of freedom) properties of some g-tori 
solutions. 



1 



1 Introduction 



In a recent paper [10] we introduced the concept of q-tori in the space of normal modes (q-space) of the 
Fermi - Pasta - Ulam model pQ described by a quadratic potential of nearby particle interactions plus a 
quartic perturbation (FPU— /3). A q-torus is an invariant torus hosting quasi-periodic trajectories whose 
limit, for /3 = 0, reduces to the product of more than one independent normal mode oscillations. In [10] the 
existence of q-tori was studied in a Cantor set in phase space, using the method of Poincare - Lindstedt 
(PL) series. Furthermore, we examined in detail the case of q-tori whose dimension is much smaller 
than the number of degrees of freedom N. Such tori were found relevant to interpreting the phenomenon 
of energy localization, which is observed in FPU trajectories resulting from particular classes of initial 
conditions. According to this, we observe that, if energy is given initially to a small set of low-frequency 
modes, the flow of energy from the initial to the remaining modes may appear 'freezed' for quite long times, 
resulting in the system exhibiting important deviations from energy equipartition. Energy localization is 
characterized by the formation of a so-called 'metastable state' |12 } 113 } 120 1 121 ( 122]. A main feature of such 
a state is the appearance of an approximate 'plateau' in the low-frequency part of the energy spectrum 
E q (a so-called 'natural packet' of modes [HI El EH]) > accompanied by an exponential tail of the energy 
distribution for the remaining modes. Important problems concern the dependence of a packet's width s 
on the specific energy e = E/N given to a system (where E is the total energy), as well as the dependence 
of the time T, up to which a metastable state persists, on e. A state-of-the-art account of numerical 
facts related to these problems is given in [19], containing also many references to earlier works. Briefly, 
numerical experiments indicate that the functions s(e) and T{e) behave differently in the a and in the f5— 
FPU models, depending also on the precise form as well as on the degree of phase coherence of the initial 
conditions characterizing packet excitation (see also the review |16j). 

Another concept related to q-tori is the concept of q-breathers, introduced in the series of works [3]-[9]. 
The q-breathers are periodic orbits, whose localization profiles share many common features with those of 
particular q-tori solutions. In fact, a q-breather can be considered as a one-dimensional q-torus. Contrary 
to q-tori, however, the existence of q-breathers is guaranteed by basic theorems on the continuation of 
periodic orbits from linear to non-linearly perturbed systems. On the other hand, the q-tori effectively 
describe a wider class of initial excitations in the FPU problem than the q-breathers. Thus, q-tori and 
q-breathers should be regarded as complementary interpretational tools of FPU localization phenomena. 
Further comparison is deferred to section [3^1 

In a paper accompanying the present one [TT], we explain in detail the numerical procedure by which 
we locate the q-tori, and we discuss the convergence properties of the corresponding PL series. In the 
present paper, on the other hand, we mainly focus on the physical aspects of our study. In particular, we 
seek to explain theoretically the energy localization properties observed in q-tori. To this end, we explore 
a concept introduced already in [10], namely the so-called propagation of 'initial excitations' in the PL 
series terms. In the context of PL series, an 'initial excitation' means a particular selection of a subset 
of modes in q-space for which we consider an oscillation with non-zero amplitude at the 0-th order of 
perturbation theory. Then, this excitation 'propagates' to new modes at subsequent orders. Similarly to 
the (5 case considered in [10] , in section [2] we prove a general proposition regarding the pattern of this 
propagation, applicable to both the a and j3 models. 

The latter proposition is exploited, in turn, in order to predict theoretically the form of energy local- 
ization profiles depending on the initial excitation. In section [3j we study such profiles and compare them 
with the profiles of FPU-trajectories arising from the energy given initially exclusively to a few modes (this 
is in contrast to the q-tori solutions in which all modes have initially some, albeit possibly quite small, 
energy). Our numerical evidence is that the localization profiles of q-tori and nearby FPU-trajectories are 
very similar. In the case of 'extensive' initial excitations (i.e. ones where the number of initially excited 
modes varies proportionally to N), we predict that the form of the energy spectrum corresponding to a 
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given localization profile depends only on the specific energy e = E/N . Furthermore, we find approxima- 
tive formulae for the form of localization profiles corresponding to various arbitrary choices of the set 

1 = 1, . . . , s of initially excited modes, based on a leading term analysis of the associated PL series. In the 
case of low-frequency packets these profiles are exponential, with a slope depending logarithmically on the 
system's parameters. We also study the case of g-breathers, where a similar law is found in [7], when the 
so-called 'seed mode' qo is taken to vary proportionally to N. All these results are verified in numerical 
computations. 

Section U] is a summary of our basic conclusions from the present study and a discussion on future 
perspectives. 

2 Analytic construction of q— tori 

2.1 The Fermi Pasta Ulam model 

The FPU Hamiltonian for a lattice of N — 1 particles reads: 



where is the k— th particle's position with respect to equilibrium and y^ its canonically conjugate mo- 
mentum. Fixed boundary conditions are defined by setting xq = xjy = 0. The cases a / 0, j3 = 0, and 
a = 0, (3 7^ are called FPU-a and FPU-/3 model respectively. 

The normal mode canonical variables (Q q ,P q ) are introduced by the linear canonical transformation 



N-1 N-1 N-1 n N-1 





k=l k=0 k=0 k=0 





(2) 



Substitution of ((2|) into ([T|) yields the Hamiltonian in the normal mode space ('g-space'): 





a 



N-l 



m 



(3) 



9=1 



q,l,m=l 




N-l 



n 



q,l,m t n=l 



with normal mode frequencies 




1 < q< N-l 



(4) 



The harmonic energy E q of each normal mode q is given by 




(5) 
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The coefficients B q i m and C q i mn are non-zero only for particular combinations of the indices q,l,m,n, 
namely 



± 

Cqimn = ^(^g±i±m±n,0 ~ 8q±l±m±n,2N) • (6) 



In the above expressions, all possible combinations of the ± signs must be taken into account. In the new 
canonical variables, the equations of motion are: 

Qq ~\~ QqQq — \/ 2N ^ ^ Bql m VL q VLlVl m QlQ m 
l,m=l 

2 jy ^ ~] Cqi mn ClqQi£l m £l n QiQ m Q n . (7) 

l,m,n=l 

2.2 Construction of q— tori by Poincare— Lindstedt series 

As a starting point for the construction of quasi-periodic solutions lying on 'g-tori', we consider first the 
trivial case a = f3 = 0. Let 

A) = {<7l, Q2, • • • , Qs}, where 1 < qi < N - 1 with qi < qj, for i < j 

be an arbitrary set of s < N modes, called hereafter 'seed modes' (in analogy to [3j, where the case s = 1 
corresponding to g-breathers is considered). The set of functions Q qi (t) = A qi cos(O gi t + 4> qi ), i = 1, . . . , s 
and Q q (t) = if q ^ 2?0i constitutes a particular solution of the linear system. The resulting trajectory 
lies on an s-dimensional torus, provided that the frequencies £l qi of Eq. ([!]) satisfy no commensurability 
relation. This turns out to be always the case if N — 1 is a prime number or log 2 N € N* [26]. If, on the 
other hand, the frequencies £l qi satisfy s' linearly independent commensurability relations (0 < s' < s), the 
trajectories lie on a 'resonant torus' of dimension s — s' , which is a sub-manifold of the original q-torus of 
dimension s. 

Passing now to the nontrivial case a ^ 0, or /3 ^ 0, we aim to define quasi-periodic trajectories lying 
again on s-dimensional tori. To this end, let iv qi , i = 1, . . . ,s be a set of frequencies with fixed values 
chosen in advance, which are incommensurable between themselves as well as with each one of the linear 
frequencies £l q of the remaining modes q ^ T>q. We then determine formal solutions Q q (t), q = 1, . . . , N — 1 
containing only trigonometric terms of the form cos(n • (cot + (/))), where n = (ni,n2, . . . ,n s ) is an s- 
dimensional integer vector and oj = (uj qi , . . . , u qs ),</>= {4> qi , ■ ■ ■ , <t>q s ) are the frequency and phase vectors 
respectively. 

According to the PL method, these solutions for FPU-q or FPU-/3 model are written as series in 
powers of a small parameter fi = a/s/2N, or fi = (3/2N respectively, namely: 

oo 

Qq(t) = Y,» k Q q k) (t), q = l,...,N-l . (8) 

k=0 

(k) 

The series terms Q q are computed step by step. The 0-th order terms are set as 

n(°)m - / A?cos(u; 9 i + <^), ifq£V 

W I 0, otherwise . { } 
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We emphasize that the amplitudes A qi , i = 1, . . . , s are unknown quantities to be specified at the end of 
the process. In the computer-algebraic program, A qi are symbols carried all along the construction of the 
PL series, while the frequencies oj q . are substituted at the beginning by their selected numerical values [TT] . 
However, according to the PL method, the frequencies uj qi must also be expressed in the form of a series 
in powers of the amplitudes A qi , namely 

oo 

"« = n ft + X> fc ai fc) 0V--->^.) • (io) 

fc=l 

(k) 

The functions ui qi (A qi , . . . , A Qa ) are polynomials of the amplitudes A qi , . . . , A qs , of order k + 1 in the 
a-case, or 2k + 1 in the /3-case. 

Substituting Eqs.(|8]) and (HOP in the equations of motion ([7]), we find the following equations to be 
solved at order k: for the FPU-a model we have 

k m 

qP+<4q?> = E E 4 m H ni ~ n2) Qt ni) 

ni=l n2=0 

N-l fc-1 

- n q E W m B qlm E Qj" ^. if^e^o 

l,m=l ^1,2—0 

n-j^ +n2~ fc — 1 

7V-1 fe-1 
Z,m=l Hi,2=0 

71 1 + "2 = A; — 1 



(11) 



while for the FPU-/? model we have 

Q q k) +^ q Q q k) = E E 4 n2 H ni ~ n2) Q ( t ni) 

ni=l ri2=0 

AT-1 fc-1 

n, E ^^n^ mn e Oi" ^ 5 ^. if^^o 

l,rn,n=l ^1,2,3—0 
iV-l fc-l 

Q q k) + n 2 q Q q k) = -n, E "i^^Q™ E Q/ (ni) ^ 2) Ql n3) , ifg££>o- 

£,m,n=l "4,2,3=0 

Til + ri9+no=fc — 1 



(12) 



By integrating either Eqs. (|lip or (|12p . secular terms of the form tsm.{ui q t) appear in the right hand side 
(when q G T>q), at even orders in the FPU-a and at all orders in the FPU—/?. The requirement to eliminate 

(k) 

all secular terms leads to an algebraic expression for the frequency correction terms co qi , i = 1, . . . , s. After 
eliminating the secular terms, direct integration of Eqs. (jlip or (|12p yields the solution of the series terms 
Qf\t) for all q = 1,...,N- 1. 

At the end of the process, we consider Eqs. (|10p as a set of algebraic equations to be solved for the 
amplitudes A qi . The existence of a solution for a fixed choice of frequency values cj qi , along with the 
convergence of the series © for that particular solution, constitutes a proof that a g-torus with the so 
chosen frequencies exists. In practice, we can hardly provide such a proof by rigorous means. Instead, we 
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work with a finite truncation of the series (jlOp . which we solve numerically to determine values for the 
amplitudes A qi , up to the approximation of the selected truncation order. Even so, a random choice of 
frequency values u) qi often results in that real solutions do not exist. In order to deal with this problem, 
we developed a numerical algorithm by which to properly choose values for the frequencies u qi , so as to 
lead to a successful determination of real amplitudes for all series terms. The same algorithm takes care 
of how to move in frequency space, in order to determine solutions corresponding to progressively higher 
and higher energy. Details on this algorithm are given in 01 

After the determination of real amplitude values A qi , we substitute these latter values in the (also 
truncated) series ([8]) . This gives finite expressions representing motions on a g-torus up to a finite accuracy. 
However, by use of some numerical indicators (like the GALI; see [23]) we provide numerical evidence on 
whether a solution found by the above approach represents indeed a good approximation to a g-torus 
solution. In fact, the GALI indicator also allows us to probe the dimension of this torus, and check 
whether the latter agrees with the dimension adopted a priori in the analytic solution. Further numerical 
tests on the convergence of the PL series are discussed in [TT] . 

2.3 Sequence of mode excitations 

It is well known that for some special choices of initial conditions, the resulting FPU-trajectories (in both 
the a and (5 models) take place on lower dimensional invariant sub-manifolds of the FPU phase space. This 
is due to the existence of discrete symmetries, which give rise to explicit low-dimensional FPU solutions. 
An extensive study on such symmetries is made in [271|28j|29] (in the latter such solutions are called 'bushes 
of normal modes'). A particular case are periodic trajectories, arising from exciting only one of the modes 
qo = iV/3, N/2 or 2N/3 in the f5 model [2] (see also [27} 129]). Solutions like the above lead, by definition, 
to energy localization, since the energy remains always distributed among a small subset of modes. 

On the other hand, as pointed out in the introduction, we find numerically that energy localization 
occurs also for sets of initial conditions not obeying any obvious or simple symmetry. We now study 
this phenomenon using the concept of 'propagation' of some initial 'excitation' in the PL series for low- 
dimensional tori. We can briefly state our main result as follows: through the study of propagation, we 
can define a hierarchy of groups of modes participating in a g-torus solution, such that all the modes in 
one group share a similar (in order of magnitude) amount of energy, while distinct groups share quite 
different amounts of energy. The hierarchy of these groups allows us to predict the whole localization 
profile via a leading order analysis of the associated PL series. However, it also allows us to characterize 
the paths of energy transfer in (/-space for FPU trajectories neighboring some g-torus solution, from an 
initial excitation up to the moment when a metastable profile is established for the FPU-trajectories. 

We now give the following definitions: 

Definition 1: A mode q is said to be excited at the n-th order of the PL scheme, iff in the series ([8]) it 
is Q q k \t) = for all k < n, and Q q n \t) / 0. In addition, Q q n \t) is called leading order term of the series (JS]). 

Definition 2: Let Dq be a set of modes excited at the 0-th order of the PL scheme according to Eqs.([9|). 
The sequence of sets V^, k = 1, 2, . . ., where consists of modes excited at the A:-th order of the PL 
scheme, is called sequence of mode excitations. 

Definition 3: Let fen be a positive constant. We call tail modes with respect to ko the modes belonging 
to the set Uk>k(P k . □ 



1 In [22], studying solutions corresponding to an initial excitation of the first normal mode, tail modes were called those 
belonging to the last third of the spectrum. This is equivalent to state that ko = [2iV/3]. 
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Definition 4: We call FPU-trajectory with initial excitation in the set of modes T>q = {gi, . . . ,q s } the 
trajectory resulting from the equations of motion ([7|) for the set of initial conditions Q q (0) = A g coscj) g , 
-fg(O) = ~~ A q Q q sin (f) q , for some set of amplitudes A q and phases cp q , if q £ Vq, and Q q (0) = P q (0) = if 

In order to determine the sequence V k produced by a particular initial excitation Vq, we first define 
the set 



$7 = {-i,+ir 



(13) 



The elements of S r are r-dimensional vectors of the form u^ r ' = (o"i, . . . , oy), where o~i = 1 or — 1, 
i = 1, . . . ,r. Furthermore, for an r-vector x = (x\, . . . ,x r ) we define ■ x as the Euclidean product 
a ( r ) . x = a\X\ + . . . + a r x r . We can now prove the following 

Proposition: Let Dq = {q±, q2, ■ ■ ■ , q s }, 1 < q\ < . . . < q s < N — 1 be the set of seed modes of an 
initial excitation yielding a formal PL solution associated with trajectories on an s-dimensional q-torus. 
Let M k be the set 



(r)Jr) I +JV- 1 



a^'q 



2N 



N- | <T (r) q 



(r)„(r) 



. a (r) € S r. ; q (r) £ p r 



(14) 



where r = r(k) with r(k) = k + 1 for the FPU-a model and r(k) = 2k + 1 for the FPU-/3 model. Then, the 
sequence of mode excitations T>k corresponding to the initial choice T>q is defined by the recursive relations 



V k = M k \ |J Vj, 



k 



1,2, 



(15) 



o<j<fc-i 



An explicit proof of the above proposition is given in [10] for the FPU-/3, and in [5] for the FPU-a. We 
note that [•] is the integer part of a number and M k represents the set of all modes for which the right hand 
side of Eqs. (jlip and (|12p is non-zero, i.e. the modes yielding some non-zero contribution to the energy 
spectrum up to the k-th order of the PL series (some modes excited at previous orders up to k might also 
belong to M k ). 

Some examples clarify the use of Eq. (|15l) : 



i) q-breathers: If we choose T>q = {qo}, the so-induced PL solution corresponds to a one-dimensional 
torus, i.e. a periodic orbit. In this case we find: 



= Uk} with q k 



rq + N - 1 
2N 



N — rqo 



(16) 



where r(k) = k + 1 in the FPU-a, or r{k) = 2k + 1 in the FPU-/?. This rule coincides with the one given 
in [9] for an arbitrary seed mode go- From Eq. (|16p we readily find that, if go is even, only even modes 



2 As rendered clear with specific examples below, in the case of g-tori one has in general Q q (0) / or P q {0) ^ for all 
q = 1, . . . , N — 1. In fact, for the modes belonging to the fc-th set T>k (see Definition 2) we have in general Q q (t) ~ 0(/j, k ) 
for all times t, including t = 0. This implies that in g-tori all modes have some energy already at t = 0. On the contrary, 
according to the Definition 4, in FPU-trajectories only the modes belonging to Do share the total energy at t = 0. In that 
sense, the use of the term 'excitation' for FPU-trajectories is literal, i.e. T>o refers to the modes excited at t = 0. The origin 
of the term 'FPU-trajectories' is that these are trajectories with initial conditions of the same type as those considered in the 
original FPU paper. 
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become excited at subsequent orders in both the a and (5 models. On the other hand, if go is odd, in the 
a model both odd and even modes become excited, while in the /3 model only odd modes become excited. 
The localization properties of solutions corresponding to (/-breather excitations will be discussed in detail 
in section 3 below. 

ii) Example with two seed modes: Suppose Vq = {51,92} = {3, 5} for N large in FPU-a. At first order 
(k = 1) it is r(l) = 2. In order to determine Mi, we consider all possible combinations of the symbols 
o-W G S 2 and g( 2 ) G V\. These are {(1, 1), (1, -1), (-1, 1), (-1, -1)} and {(3, 3), (3, 5), (5, 3), (5, 5)} re- 
spectively. Then from Eq. (fT41) we find M\ = {2,3,5,6,8,10}. Since Vq = {3,5}, from Eq. (fT5l) we have 
V 1 = Mi\V Q = {2, 6, 8, 10}. Repeating the above procedure for k = 2, we find M 2 = {1, 3, 5, 7, 9, 11, 13, 15} 
and T>2 = M2 \Vq VJV\ = {1, 7, 9, 11, 13, 15}. In the same way we proceed to subsequent orders k = 3, 4, . . . 
In the FPU-/? one follows the same steps, but for r{k) = 2k + 1. 

iii) Excitation representing a low-frequency packet of seed modes: Let us consider as an initial excitation 
the packet of modes Vq = {1,2,3,4} for ./V large. Following the same procedure as above, in the FPU-a 
we find the sequence of excitations V\ = {5, 6, 7, 8}, V2 = {9, 10, 11, 12}, etc. We notice that at each order 
modes are excited in groups. In the same way, the j3 model yields V\ = {5, . . . , 12}, T>2 = {13, ... , 20}, 
etc. 

These propagation rules can be generalized for s-dimensional g-tori corresponding to low-frequency 
packets of modes. The seed mode excitation Vq = {1, 2, . . . , s} generates V k = {ks + 1, . . . , (k + l)s} in the 
FPU-a, and V k = {(2k — l)s + 1, . . . , (2k + l)s} in the FPU-/3, with k > 1. Furthermore, we observe that 
V k = I^fc-i U ^2k- As shown in the next subsection, these rules imply that the resulting g-tori solutions 
exhibit exponential energy localization profiles. 

iv) Excitation representing a high-frequency packet of seed modes: As an example, let us consider 
V = {28,29,30,31} in the TV = 32 dimensional chain. In the FPU-a we find V x = {1,...,8}, then 
T>2 = {20, . . . , 27}, etc. The general rule, by setting the last s modes Vq = {N — s, . . . , N — 1} as seed 
modes, is: V k = {(k - l)s + 1, . . . , (k + l)s}, if k = 2n+l and V k = {N-(k + l)s, . . . ,N- (k- l)s- 1}, if 
k = 2n. By the same way, in the FPU-/3 we find T> k = {N — (k + l)s, . . . , N — (k — l)s — 1}, V/c. Comparing 
the two models, we see that V{ = V% k . 

v) Discrete symmetry solutions: Suppose Vq = {A^/2}, or Vq = {2iV/3}. We then find V k = Vq 
for all k = 1,2, .. . in both models, while in FPU-/3 the condition V k = Vq holds also for Vq = {A r /3}. 
The g-breather solutions for these cases correspond to the 'nonlinear normal modes' of the FPU system 
[2"1 1271 [29], that coincide with the FPU-trajectories resulting from the same seed mode. 

We note finally, that for g-breathers we have a general relation connecting the sequences of excitations 
V k in the a and in the (5 model, starting from the same seed mode. Namely, from Eq. (|16p we find that 
V^ = V^ k . In words, the mode excited at the &-th order in the FPU-/3 is the same as the mode excited at 
the 2fe-th order in the FPU-a. This relation holds also for excitations of small packets around A^/4, A^/2 
and 3N/4, but it does not hold in the case (iii) (low- frequency packets of modes). In fact, for an arbitrary 
excitation we have M% k = ikff , but we only have £>f = V% k provided that M^,_ 1 n M% k = 0. 

2.4 Leading order terms 

We close this section by giving an expression for the leading order terms of the solution Q q (t), q £ V k (see 
Definition 1). This will be used in section [3] in order to predict the forms of energy localization profiles of 
g-tori. 
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Let Vq = {qi, Q2, . . . , q s } be a 0-th order excitation set, G Vq an r-vector in V r and w„ = 
{uj ni , . . . ,u nr }, <pn^ = (4>m, ■ ■ ■ , 4>n r ) its associated frequency and phase vectors respectively. Some basic 
features of the PL construction arising in accordance with the Definitions 1 - 4 of subsection 12.31 are: 

i) If q £ X>fc, then Q q k \t) is a leading order term of the series (jHJ). 

ii) The only modes that admit frequency corrections are those in Vq, for the rest holds uj q = Q q , q ^ Vq. 

iii) At the fc-th order, the expressions for Q q k \t) contain divisors which are the products of k factors 
of the form — (crW" 1 )) < k. Since all frequencies are fixed in advance to incommensurable 
values, such divisors, albeit possibly quite small, can never become equal to zero exactly. 

We determine the solutions of Eqs. ljlip and (|12p . for the leading order terms Q q k \t), starting by the 

0-th order solution Q ( q \t) of Eq. ©. For convenience, in subsequent formulae we use exponential rather 
than trigonometric expressions. In both the a and /3 models, for q G V^ we find (seeUJ 

Q(*>(t)= y ^(nWjJCf)^))^^) (17) 



where r(k) = k + 1 for the FPU-a and r(k) = 2k + 1 for the FPU-/3. In the above expression, the factor 

(k) 

lZ q is given by 

and idq^ by the recursive relation 

4S(- (r) )= E E «(« (r(k)) )«(^ (19) 

il,2=0 rn i &T> li 

l 1 +l 2 = k-l i = l,2 

with 

n (r(fc)) = (ni,...,n r(/l) ,n r(il)+1 ,...,n r(fc )) 



in the a case and 

k-1 



^M T) )= E E ^)(n(KM) )£ (^) (n (K^)) )£ g (n 



h, 2,3=0 rriieVi. 
i l + ! 2+ i 3= fe - 1 i=l,2,3 



with 



(r(fc)) 



(ni, . . . , "-r(h)' n r(/i)+l) ■ ■ ■ ) n r(/i)+r(« 2 )' n r{h)+r{h)+li ■■■i n r(k)) 

" v ' V «, ' V „ ' 



in the /3 case, setting, in both cases, K, q = 1 at k = 0. The terms £m, m 6 Pi entering the expressions 
(Trail and flU]) are 

£(0 (n (K0)) = (21) 
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for k > 0, or £ q = 1 at k = 0. Both quantities IC q -a, and K> q .n are polynomials of degree k — 1 in the 

terms £m , Z = 0, 1, . . . and can be computed iteratively. Explicit expressions for the first few orders of the 
mapping (|19p are given in[D 

3 Numerical examples. Localization profiles 

In this section we provide various numerical tests on the dynamics of the g-tori and their neighboring 
FPU-trajectories. In these, we compare computer-algebraic series calculations with purely numerical 
calculations. In particular, we examine and correlate the energy spectra in both cases. Furthermore, we 
derive the form of the localization spectra for particular g-tori solutions and check up to what extent 
various features predicted by theory appear in the localization spectra of the numerically integrated orbits. 
In H] we give a table with the frequencies and amplitudes used throughout all the examples below, in order 
to allow for reproducibility of the results. 

3.1 A q— torus example 

As explained in subsection [221 hi numerical computations we work with finite truncations of both series ([S]) 
and ()10p . Furthermore, for a specific choice of frequencies ui qi we solve numerically the truncated Eqs. (|10p 
in order to specify the corresponding numerical values of the amplitudes A qi . Then, we substitute A qi into 
the series (|SJ) truncated at the order ko. 
In the sequel, we denote by 

Qq L,k °(t) = QfHt) + /uQ? } (*) + • • • + v^Qt °\t) 

q = 1, . . . , N — 1 the approximate solutions of Q q (t) obtained by the above procedure. By differentiating 
with respect to time, we also specify a truncated series for the canonical momenta P q ,k °(t) = Q q L,k °{t). 
Having both expressions at hand, we specify initial conditions Qq L ' °(0) and P^ ,k °(0) which, in this 
approximation, should lead to motion on a torus with the prescribed frequencies. To verify the above 
statement, we integrate numerically the orbital resulting from such initial conditions and obtain Q q Ln,k ° (t), 
P q PLn ' °(t), in order to make various comparison tests between the PL and their numerical realization (PLn). 
A main comparison refers to the normalized averaged energy spectra in a time interval 2T/3 <t<T. We 
define the average harmonic energy of the q-th mode as 

E q (T) = | f T E q (t)dt (22) 

1 J2T/3 

where E q is given by Eq.([5|). By introducing the re-scaled wavenumber k = q/N, we plot the normalized 
averaged energy spectrum e K = E q /E versus k, where E is the total energy given to the system. Such 
plots allow us to compare the averaged energy spectra eK L,k °, eK Ln ' k °, and also the spectra e K of FPU- 
trajectories, for various solutions corresponding to N fixed, but also iV varied, which is a case of interest 
in the solutions arising from so-called 'extensive' sets of initial conditions (see below). 

As a first example, we consider a 4-torus solution in the FPU — a model with N = 32, a = 0.33, when 
the initial excitation refers to the first four lowest frequency modes, i.e. T>o = {1,2,3,4}, with the chosen 
frequencies given in the first group of the table in HI The truncation order is ko = 1 1 and all phases are 
set as 4> qi =0. Solving numerically Eqs. lfTU]) we specify the amplitudes A±, A2, A3, A4 (seed]), which result 
to E = 0.000182466. 

3 by a 4-th order symplectic kinetic-potential splitting algorithm 
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Figure 1: (a) Normalized averaged exponential spectra e K versus k = q/N for the system with T>o = 
{1,2,3,4}, a = 0.33, E = 0.000182466, N = 32, for a total time T = 10 6 . The blue spheres correspond to 
a g-torus construction with the truncated PL series, denoted Qq L,11 (t). The orange triangles correspond to 
a numerical solution, denoted Qq Ln ' ll (t), obtained by numerically integrating the FPU equations of motion 
with initial conditions Q q ' n (0), P q ' (0). The continuous line is log(e K ) = 3.49545 — 36.6613k — 2 log (re), 
(b) GALI indices (see text) for k = 2, 3,4, 5 calculated for the orbit of the orange spectrum in (a). 

FigQTa) shows a comparison of the normalized averaged spectra e^ L (blue) and e^ Ln (orange) at 
integration time T = 10 6 . The continuous line is obtained by the method described in subsection 13.21 
below, which yields a prediction for the mean exponential slope of the energy localization profile based 
on the leading order terms of Eq. (]17|) for an initial excitation corresponding to the first s lower frequency 
modes. 

Since we only compute the numerical coefficients of the truncated PL series solution, we have no 
rigorous demonstration of the convergence of the series. In [TT] we will discuss several criteria examining 
this question. Here, instead, we rely, as in |10| . in the use of the Generalized Alignment Index (GALI) 
method |23 |. 124 1 [25], which is a numerical method allowing us to determine the dimension of a stable torus 
by the temporal behavior of a set of indicators (G2, G3, . . . ) computed via an integration of the variational 
equations of motion along with the original ones. FigQJb) shows in double logarithmic scale the evolution 
of GALI indices for k = 2, 3, 4, 5. The indices for k = 2, 3 converge very fast to a constant value, while 
G4 oscillates around its average value after an initial decay lasting for a rather long time, t = 3 • 10 7 . On 
the other hand, at times greater than the time of stabilization of G4, the index G5 clearly decays as 1/t, 
up to at least the time 10 9 . Following |24^ 125] . we conclude that the numerical orbit Qq Ln,ll (t), Pq Ln ' ll {t) 
lies on a torus of dimension s = 4. In fact, since we only have a finite precision in the initial conditions, a 
more precise statement, is that the orbit follows a g-torus dynamics for times up to 10 9 , so that no large 
scale chaotic diffusion phenomena can occur for this orbit within this timescale. 

3.1.1 Comparison with an FPU trajectory. Stages of dynamics 

FPU-trajectories with initial conditions close to (/-tori appear to exhibit, on the average, similar localization 
phenomena as the g-tori for quite long integration times. However, as in [22], we also find that the FPU- 
trajectories undergo some well distinct 'stages of dynamics' that further characterize the time evolution of 
their average normalized energy spectra. 

As an example, we consider the FPU-trajectory with initial conditions Q q (0) = A q , for q = 1,2,3,4, 
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Figure 2: (a) Time evolution of the normalized averaged energies E q (t)/E, in double logarithmic scale, for 
the FPU-trajectory arising from the excitation of the first s = 4 modes with initial conditions Q q (0) = A q 
if q = 1,2,3,4 and Q g (0) = else, and system's parameters as in Fig (TJ for the (/-torus solution, (b) 
Comparison of the normalized averaged energy spectra e K versus k of the (/-torus solution of Fig(TJ (blue) 
and the corresponding FPU-trajectory (orange) over a total integration time T = 10 6 . 

and Q g (0) = 0, for q = 5, . . . 31, P q (0) = 0, q = 1, . . . , 31 in the FPU-a system with the same parameters 
of the (/-torus solution of FigUJ The initial distance in phase space between the FPU and the g-torus 
trajectories is of order 0(/x), where in this case fi ~ 4 x 10~ 2 . 

Fig|2^a) shows the evolution of the normalized averaged energies E q (T)/E for all the modes q = 
1, . . . , N — 1 of the above FPU-trajectory. We observe that the modes q = 1, 2, 3, 4, which shared exclu- 
sively the total energy at T = 0, continue to share most of the energy at all subsequent times T. As a 
result, E q {T)/E remains practically constant for these 4 modes, exhibiting only small variations (of order 
10~ 4 ) which correspond to the energy gradually transferred to the remaining modes. Such variations are 
unobservable on the scale of Figj2] in the corresponding curves for the modes q = 1, 2, 3, 4. 

However, the energy transfer is observable by plotting the time evolution of E q (T)/E for all other 
modes. In fact, we observe that this transfer takes place in a rather short time interval, of order T = 10 3 
for the modes q = 5, 6, 7, 8, and T = 10 2 for the remaining modes. Following the terminology introduced 
in |22j . this phenomenon corresponds to the so-called first stage of dynamics in the FPU-a system. In 
[22] it was stressed that the transfer of energy during this stage appears to be a phenomenon explained by 
regular ('integrable-like') dynamics, since a similar phenomenon is observed in the Toda lattice, which is 
an integrable system. Furthermore, the energy growth during the first stage of dynamics is characterized 
by a power law growth with respect to time, a process called 'secular avalanche'. 

It is important to point out that, also here in the case of FPU-trajectories, the modes form well 
distinguished groups. These groups correspond exactly to the groups T>k of the sequence of excitations 
generated by the (/-torus solution of FigfTJwith initial excitation T>q = {1,2,3,4}, namely we have: T>\ = 
{5, 6, 7, 8}, T>2 = {9, 10, 11, 12}, or, in general, = {ks + 1, . . . , (k + l)s}. It is observed that, during 
this secular avalanche process, for all modes q ^ T>q in the same group T>k the quantity E q /E grows in time 
by almost the same power-law, i.e. the energy spectra behave as E q (t)/E ~ t Cq ( k \ where the exponent is 
almost constant within a group: c q (k) = c(k) + e q , with \e q \ « 1 for all q £ T>k- 

On the other hand, at the end of the first stage of dynamics the normalized averaged energy spectrum 
for the FPU-trajectory appears to be stabilized to a form persisting at least up to the end of the numerical 
integration, at T = 10 6 . This stabilized spectrum displays an exponential profile, which is strikingly 
similar to the (/-torus profile shown in FigQJa), as can be seen by a direct superposition of the two spectra 
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(Figj2]||b)). In particular, we observe that the modes in the spectrum of the FPU-trajectory, during this 
metastable state, remain separated in the same groups, exactly as in the g-torus solution. 

We emphasize that this effect cannot be considered as a trivial consequence of the proposition of 
subsection 12. 3\ since the latter only concerns a formal aspect of the PL series construction for g-tori, while 
the separation of the modes in the case of the FPU-trajectory is a dynamical phenomenon, characterizing 
the energy transfer process. The fact that, the final groups defined in either case are the same, indicates 
some deeper connection between g-torus and FPU-trajectory solutions. However, although the FPU- 
trajectories appear to follow effectively a g-torus dynamics up to a long time, contrary to g-tori, most 
FPU-trajectories enter eventually into the stage of approach to energy equipartition (the so-called 'second 
stage of dynamics' in [22]). 

3.2 q— Tori low frequency packet solutions and exponential energy localization 

The g-torus as well as the FPU-trajectory examined in the previous subsection are examples of a class of 
solutions of particular interest in the literature (see [TBI \T7\ [TBI f!9] ) , namely solutions corresponding to 
the initial excitation of a packet of low-frequency modes. In fact, it is numerically found that, for values 
of the specific energy beyond some threshold, low-frequency packets of modes are formed naturally, even 
if initially is excited only one, e.g. the q = 1 mode. Some questions of central interest in the literature 
concern the dependence of: i) the width of natural packets ii) the exponential slope of the energy spectrum 
of the remaining modes, on system's parameters E, N etc (see [16] for a review). 

In the sequel we examine the form of localization profiles for g-torus solutions associated with an initial 
excitation of a low-frequency packet of modes. We obtain theoretical results based on a leading order term 
analysis of the PL series for g-tori of subsection 12. 41 Furthermore, we compare these results with ones 
found numerically for FPU-trajectories with a similar initial excitation. 

Our main result for the FPU-a can be stated as follows: in subsection 12.31 it was mentioned that, for 
g-tori, an initial excitation T>q = {1, 2, . . . , s} in the PL series leads to the sequence of mode excitations 
T^k — {ks + 1, . . . , (k + l)s}- Starting, now, from the median mode in the group T>k, i-e. the mode 
qmid = ks + [s/2], we can obtain estimates of the size of the leading term Qq^ id given by Eq. (fT7|) and derive 
estimates on the magnitude of the harmonic energy Eq^ id , which is hereafter denoted by E^> . Then, we 



where M = s/N is the fraction of initially excited modes with respect to the total number of modes. The 
derivation of Eq. (|23|) is given in [H 
Normalizing Eq. (|23|) we obtain 



normalized energy profiles of q-tori remain unaltered as N increases. 

This prediction becomes hardly possible to test by a direct construction of the g-tori solutions via 
PL series, because as ./V increases, we quite soon encounter the limits of computer memory required for 
storing the coefficients produced by the computer-algebraic program. However, taking into account the 
evidence presented in subsection I3.1.1[ that FPU-trajectories with the same initial excitations as g-tori 
exhibit similar localization profiles, we can test numerically the extent up to which the invariance of the 
normalized averaged energy spectrum holds for FPU-trajectories as well. 

Such a test is made in Figj3j In (a) we give the normalized averaged energy spectrum e K as a function of 
the re-scaled wavenumber k = q/N for an FPU-trajectory with a = 0.33, M = 1/8, and e = 1.5625 x 10~ 8 . 



have: 




(23) 



log e K ~ -j^n -\- 2 log k log(A 1 / 2 M 3 iV) (24) 
where A = a 2 e/vr 4 M 4 . The main prediction is that if a, e, and the fraction M = s/N are kept fixed, the 
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Figure 3: Normalized averaged energy spectra e K versus k for FPU-trajectories (orange spheres) in both 
the FPU-a and FPU-/3 models, keeping the specific energy g and the ratio of initially excited modes M 
fixed. The parameters are: (a) FPU-a model with a = 0.33, M = 1/8, e = 1.5625 x 1CT 8 . (b) FPU-/3 
model with /3 = 0.3, M = 1/16 and e = 1.5625 x 10 -6 . In (a), the blue line in all panels corresponds to 
the theoretical prediction of Eq. (j24p based on the leading order term analysis for (/-tori. This is given by 
loge K = —56.27k + 4.46 + 21ogK. In (b) the blue line corresponds to a similar prediction for the FPU-/3, 
given by Eq.®- We find loge K = -78.64k - 1.204. 

We computed these trajectories by progressively increasing N, namely N = 64, 128, 256 and 512. The 
energy spectra of Figi3]are all evaluated at T = 10 6 . The solid line corresponds to the fitting law of Eq. 
(|24p . which, for the adopted parameters, takes the form indicated in the figure caption. The main remark 
is that the same line fits all re-scaled normalized spectra, for different N (while the fraction of initially 
excited modes M = 1/8 as well as the specific energy e = 1.5625 x 10 -8 are kept constant). 

A question of central interest regards the upper limit in the specific energy for which the normalized 
spectra for the FPU-trajectories continue to exhibit exponential localization. Eq. (j23p allows us to obtain un 
upper limit, by requiring that A = a 2 e/-7r 4 M 4 < 1. However, as we approach the upper limit e = 7r 4 M 4 /a 2 , 
the analysis based on only the leading order terms of the PL series ceases to be valid, since important 
contributions to the energy spectrum are made also by the higher order terms in each mode's series 
expansion of Eq. (jHJ) . 

The condition A < 1 implies s > / 4 iV. Thus Eq.(|23p applies when the initially excited packet 

has a width larger than the width of the so-called natural packets |14[ \T7\ \TE[ \19\ . In the case of natural 
packets, it is found that by the excitation of a number of low frequency modes satisfying s < a 1 / 2 e 1 / 4 A^, 
the so resulting energy spectrum exhibits a plateau of width a 1 / 2 e 1 / i N (larger than the initially excited 
modes). Furthermore, there is evidence that the slope a of the exponential energy localization profile 
e K ~ exp ( — o • k) depends linearly on a~ 1 / 2 e -1 / 4 [20]. This is in contrast to the slope which refers to g— tori 
solutions of Eq. (|2~3"|) . that depends logarithmically on [a~ 1 / 2 e~ 1/ ' 4 ]~ 4 and points out that different choices 
in the fraction of the initially excited low-frequency modes result in different exponential laws. 

g—Tori for low— frequency packets in the FPU— /3 For the sake of completeness we report the results 
of our previous work [10], concerning exponential energy localization in the FPU-/3 model. Assuming the 
initial excitation to be T>q = {1,2,... , s}, the sequence of mode excitations in the PL series at the orders 
k = 1, 2, ... is V k = {(2k - l)s + 1, . . . , (2k + l)s}. Estimating the size of the median mode q m id k 

= 2ks 

in each group T>^ via Eq. (|17j) . we are lead to an estimate for the energy spectra of g-tori with the above 
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excitation, namely 

Em = li < 25 > 

where M = s/N . Using re-scaled variables as in the a case, Eq. (|25|) takes the form 

l0geK ~ ~M~ K ~ l ° gS ^ 

where A = /3e/-7r 2 M 2 . We find a similar result as in the a case, namely Eq. (|26p implies that by keeping 
both the specific energy e and the fraction of excited modes M fixed, while N increases, the normalized 
exponential profile remains invariant. Again, in order that the analysis be valid, one must have s > 
-/V(/3e) 1//2 , implying that the initial excitation should be in a regime quite different from that of natural 
packets. In fact, in the present case as well, Eq. (|26p describes correctly the localization profile provided 
that A << 1 (see also (10]). 

As a numerical test of the above predictions we use again numerical computations based on FPU- 
trajectories rather than exact q— tori solutions. The four panels of Fig 12(b) show the normalized averaged 
spectra e K of FPU-trajectories, along with the predictions of Eq. (|26|) . for the fixed values M = s/N = 1/16, 
and e = 1.5625 • 10 -6 . We observe again that the spectrum remains practically invariant with increasing 
N. 



3.3 Localization patterns for arbitrary initial excitations 

So far, we focused on g-tori, and their neighboring FPU-trajectories corresponding to initial excitations in 
the low-frequency part of the spectrum. However, it is possible to see that energy localization appears also 
in cases where the initial excitation has quite different features than in the case of packets of low-frequency 
modes. In particular, we will examine now q— tori solutions in the FPU-a system produced by an initial 
excitation Dq consisting of a small set of s modes arbitrarily distributed in g-space. We give several such 
examples, in which we vary s, N, E, as well as T>q. As in the previous example of Fig{TJ in all present cases 
we compare the normalized averaged energy spectra e^ L obtained with the PL series, with the ones e^ Ln 
obtained by numerical integration of the equations of motion for the initial conditions Qq (0), Pq L {0), 
q = l,...,N-l. 

Example 1: evenly distributed initial excitation. In the FPU-a with N = 32, we construct a 
g-torus starting from the 0-th order excitation Dq = {1, 11, 21, 31}, when the frequency values for oj\, wu, 
^21, <^3i are chosen as in the second group of|H The truncation order here is ko = 11. Solving numerically 
Eqs. fflOl) we specify the values of the amplitudes A\, An, An and Aj,\. The total energy is E = 0.001563. 

Figl^a) shows the normalized averaged energy spectrum for the above g-torus solution. Energy lo- 
calization is manifestly present, since we observe the formation of four peaks of the energy spectrum 
around the seed modes 1,11,21 and 31. The localization pattern is readily understood by computing 
the sequence of mode excitations deduced by the proposition of subsection 12.31 Namely, we find 
£>i = {2,10,12,20,22,30}, V 2 = {3,9,13,19,23,29}, V 3 = {4,8,14,18,24,28}, £> 4 = {5,7,15,17,25,27} 
etc. We observe that consecutive modes, adjacent (on either side) to the initially excited ones, are excited 
at subsequent orders of perturbation theory. Thus, starting for example from the mode q = 11, the modes 
q = 10, 12 are excited at first order, q = 9, 13 at second order, etc. This explains the formation of the 
peaks in the spectrum. In fact, the pairs {10, 12}, {9, 13}, etc. share quite similar energies, corresponding 
to excitation amplitudes 0(h), 0(fi 2 ), ... As a result, the local form of the energy spectrum on either side 
of one peak is exponential. 

Finally, as evident in FigJ^a), we find a very precise agreement between the normalized spectrum 
corresponding to the analytical solution E^ L , and the one Eg Ln obtained by numerical integration of the 
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Figure 4: Normalized averaged energy spectra E q versus q, deduced by the PL series g-torus solution (blue 
spheres), or by a numerical integration (PLn) with initial conditions on a g-torus (orange triangles, see 
text), (a) N = 32, a = 0.33, V = {1,11,21,31}, ( b) N = 64, a = 1 with V = {60,61,62,63}. (c) 
N = 128, a = 1 with V = {63,64,65}. (d) N = 128, a = 1 and V = {94,95,96,97,98}. The chosen 
frequency values, and resulting amplitudes and energies in each case are given in|3J 
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initial conditions on the (/-torus. This fact indicates that at the truncation order ko = 11 the solution has 
converged to a good accuracy. 

Example 2: initial excitation in the high frequency part of the spectrum. We consider a 
(/-torus solution found by PL series in the case N = 64, a = 1, when Vq = {60, 61, 62, 63}, while the choice 
in the frequencies and the resulting amplitudes are shown in the third group of [U The truncation order of 
the series is ko = 8 and the energy is E = 0.000883. 

Figl^b) shows the resulting normalized averaged energy spectrum for the above (/-torus solution. This 
displays several features similar to the case of a low-frequency excitation. Namely, we observe the formation 
of groups of consecutive modes sharing a similar amount of energy. The sequence of mode excitations in 
this case turns out to be V x = {1, . . . , 8}, V 2 = {52, . . . , 59}, V 3 = {9, . . . , 16}, V A = {44, . . . , 51}, etc. 

Finally, we note again the exponential fall of the energy along two separate branches of the spectrum, 
namely a low-frequency and a high-frequency branch. 

Example 3: excitation in the middle part of the spectrum. We consider the middle modes initial 
excitation Dq = {63,64,65} in the FPU-a with a = 1, N = 128, and with frequencies and amplitudes 
displayed in the fourth group of [H The PL series are truncated at ko = 16 and the energy of the system 
is E = 7.63469 x 10~ 5 . 

In the corresponding normalized averaged energy spectra Eg L /E and E^ Ln /E (FigJ^c)) we observe 
that three energy peaks are formed: by the lowest, the highest and the middle modes of the spec- 
trum. We find that the sequence of mode excitations here is V\ = {1,2,126,127}, T>2 = {61,62,66,67}, 
£> 3 = {3, 4, 124, 125}, £> 4 = {59, 60, 68, 69}, etc. 

Example 4: excitation in the 3/4 part of the spectrum. We consider, as before, N = 128, 
a = 1, and an initial excitation Dq = {94,95,96, 97,98}. The chosen frequency values and the so- 
resulting amplitudes are shown in the fifth group of HI The truncation order is ko = 9 and the energy is 
E = 0.000649478. 

At subsequent orders, we now find the sequence of mode excitations T>\ = {1,2,3,4} U {60, . . . ,68}, 
£> 2 = {26, . . . , 38} U {90, . . . , 93} U {99, . . . , 102}, V 3 = {5, . . . , 8} U {56, . . . , 72} U {120, . . . , 127}, £> 4 = 
{22, . . . , 25} U {39, . . . , 42} U{86, . . . , 89} U [103, . . . , 106}, etc. This leads to the localization pattern shown 
in Figged). 

Generalization: We now generalize results on the localization patterns formed by initial excitations 
of packets, around the locations in (/-space corresponding the one fourth, half and three thirds of the 
spectrum. We suppose that each packet is of the form Dq = [kq — e, kq + e], having a width equal to 
s/N = 2e, e << 1 in the normalized (/-space k = q/N € [0, 1]. In order to specify the sequence of mode 
excitations at subsequent orders, we first specify the sets defined in Eq. (fl4j) . whereby the sets 
are immediately derived by the relation = Mfc\Uo<j<fc-iMj. Examining in detail the case of excitations 
around the re-scaled wavenumbers kq = 1/4, 1/2, or 3/4 we have: 

1) For kq = 1/4 we find = [0, e(k + 1)] U [1/2 - e(k + 1), 1/2 + e(k + 1)] U [1 - e(k + 1), 1], if 

k = 2n + 1, or = [1/4 - e{k + 1), 1/4 + e(k + 1)] U [3/4 - e(k + 1), 3/4 + e{k + 1)], if k = 2n. The 

resulting localization pattern displays three peaks around k = 0, 1/2, and 1, produced at odd orders, and 
two peaks around 1/4 and 3/4, produced at even orders. The total pattern is shown in Figj5]^a). In this 
figure, the ordinate in all panels indicates the order k of the PL series at which the corresponding mode, 
of wavenumber k = q/N, is first excited. In fact, according to the leading order term analysis of the PL 
series discussed above we have log e K ~ k. Thus, the patterns shown in all panels of Figj5]are similar to the 
normalized averaged energy spectra for the corresponding excitations when plotted in semi-logarithmic 
scale. 
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Figure 5: Propagation of modes for an initial excitation Dq = [kq — e, kq + e] when (a) kq = 1/4, (b) 
Ko = 1/2 and (c) kq = 3/4. In all panels, the ordinate yields the order k of perturbation theory at which 
the mode q is first excited. Since in the exponential energy localization regime we always have log e K ^ k, 
the patterns shown in the present panels are similar to the energy spectra for the same excitations plotted 
in semi-logarithmic scale. 

2) For K = 1/2 we find AfW = [0, e(k + 1)] U [1 - e(k + 1), 1], if A; = 2n + 1, or M<*) = [1/2 - e(k + 
1), 1/2 + e(k + 1)], if k = 2n. Thus, we have two localization peaks around k = and k = 1 created by 
contributions at odd orders of the PL series, and one more peak around k = 1/2 for contributions at even 
orders. The overall localization pattern is shown in Figj5jb). 

3) For kq = 3/4 we find similar results as in case (1), since AfW = [0, e(k + 1)] U [1/2 - e(k + 1), 1/2 + 
e(k + 1)] U [1 - e[k + 1), 1], if k = 2n + 1 and k > 1, and = [1/4 - e(k + 1), 1/4 + e(k + 1)] U [3/4 - 
e(k + 1), 3/4 + e(k + 1)], if k = 2n. In fact, the only difference with respect to case (1) concerns the order 
k = 1, for which it turns out that = [0, 2e] U [1/2 — 2e, 1/2 + 2e\. The resulting localization pattern is 
shown in FiglSJc). 

3.4 g-breathers and FPU— trajectories 

The existence, stability, and energy localization properties of g-breathers have been studied extensively 
in [3]-[9]. For low energies, it is found that FPU-trajectories with an initial excitation of the same 'seed 
mode' and energy as stable g-breathers exhibit similar energy localization phenomena as the g-breathers 
for quite long integration times. 

In the sequel we examine g-breathers as a particular case of one-dimensional g-tori constructed via 
PL series. At any fixed order k, the number of terms in the resulting series is substantially smaller for 
g-breathers than for g-tori of any other dimension s > 1. This fact allows us to construct the series up 
to a very high order (in the case N = 32 we were able to compute examples of PL series for g-breathers 
up to the truncation order ko = 250). Even so, for reasons explained in the accompanying paper |llj . the 
convergence of the resulting series is very slow, and in practice we obtain little gain in precision after a 
truncation order near ko = 50. In most of our tried examples, the precision achieved for the computation 
of initial conditions on a g-breather using PL series is typically 4 to 6 significant digits, while in some cases 
we reach 10 (see below). However, we also use these numbers as an initial guess in order to determine 
more digits in the initial conditions via a root-finding technique. This is possible especially for g-breathers, 
which are periodic trajectories, while we cannot use such technique in the case of g-tori of dimension higher 
than one. Using the Newton-Raphson method, we specify initial conditions on g-breathers up to 12-14 
significant digits. 

As a first example, we consider the classical case where we excite initially the mode go = 1 in the 
FPU-a. The sequence of mode excitations is g^ = k + 1. Fig J6] refers to a calculation for N = 32, a = 0.33 
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Figure 6: FPU-a system with a = 0.33, N = 32 and E = 0.00016635. In Panels (a), (b) and (c) is shown 
the evolution of the normalized instantaneous spectra E q (t)/E of modes q = 1, . . . , 8 for the PL solution, 
the numerical integration of the initial condition Q^ L (0), P^ L (0) and for the FPU-trajectory respectively 
(Ei/E has the highest value and the rest progressively decrease). In panel (d) is the exponential profile 
of the normalized averaged energy spectra e K versus k of panel (a) with blue spheres, of (b) with orange 
triangles, of (c) with green triangles, as well as the ^-breather found by Newton-Raphson (NR) with black 
spheres. The continuous line is the fitting law (|27p . 
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and E = 0.00016635, for the truncation order ko = 51. The choice of u±, as well as the resulting amplitude 
A\, are given in[U 

FigJ5fa),(b) and (c) shows the evolution of the normalized harmonic energies E„(t)/E, for q = 1, . . . , 8 in 
three different computations. Namely, in (a) we compute E q (t) by the analytical solution Q q L ' 51 {t) as found 
by the truncated PL series. In (b), we integrate numerically the initial conditions Q q ,51 (0), P q ,51 (0). 
Finally, in (c) we consider an FPU-trajectory rising by the initial condition x n (0) = A\ sm(irn/N), y n (0) = 
0, n = 1,...,N- 1. 

The main remark, by a comparison of the three panels, is an important difference in the temporal 
behavior of the (/-breather from that of the corresponding FPU-trajectory. Namely, the energies E q (t), 
q = 1, . . . , 8 remain practically constant in the case of the (/-breather, while they behave as quasi-periodic 
functions in the case of the FPU-trajectory. In fact, the energies E q (t) for the latter oscillate around 
mean values following closely the energy values of the (/-breather solution, but with an amplitude causing 
variations of more than one orders of magnitude. 

We note that no considerable difference exists in panels (a) and (b), i.e. between the time evolution 
of the energies E q computed via the analytical (/-breather series solution and its numerical realization. 
In fact, we find that this trajectory returns to its initial conditions after the time T = 2ir/ui, with uj\ 
specified by the PL method, up to 10 significant digits. This serves as an estimate of the overall precision 
in the analytic determination of the periodic orbit using PL series at the truncation order ko = 51. 

An even greater precision can be reached, however, using the values Q q L ' 51 (0), P q L ^ l {Q) as initial 
guess values for a numerical (Newton-Raphson) determination of the periodic orbit. FigH^d) shows a 
comparison of the averaged normalized energy spectra in all four computations, namely (i) Q q L ' 51 (t), (ii) 
Qq Ln ' 51 (t), (iii) the FPU-trajectory, and (iv) the periodic orbit with initial conditions as determined by 
the Newton-Raphson. The fitting line in Figj6]^d) corresponds to the law 

E q = 1 q - l q 2 E l , where 7 = a 2 JV 4 e/"" 4 > (27) 

given in [3]. We conclude that, in agreement with [22], after a short transient time, the FPU-trajectory in 
this case acquires an energy spectrum similar to the one of the corresponding (/-breather solution. 

FiglT] refers, now, to a different (/-breather example, in which we choose to initially excite a mode 
at a randomly chosen position in (/-space, namely go = 25, in the system with a = 0.33, N = 32 and 
E = 0.00465079. Again W25 and A25 are found in the table of 131 while the truncation order here is ko = 60. 

The sequence of mode excitations derived from Eq. (|16p up to the 10-th order of perturbation theory 
is q 1 = 14, q 2 = 11, q 3 = 28, g 4 = 3, q 5 = 22, q 6 = 17, q 7 = 8, q 8 = 31, q 9 = 6, and q w = 19. 

FigfT^a) shows the time evolution of the normalized energies E q (t)/E of the first seven modes in 
the sequence of excitations, namely qo, . . . ,q$, for the PL solutions Q q L,m (t) (blue) and their numerical 
realization Q q Ln,G0 (t) (orange). Differences observed between the two solutions gradually increase (on a 
logarithmic scale) as we go to higher-order modes in the excitation sequence. 

FiglT^b) now shows the normalized averaged spectra for the solutions (i) Q q L,6 °(t), (ii) Q q Ln ' 60 (t), 
and (iii) an FPU-trajectory rising by the seed mode excitation qo = 25. We clearly see again that the 
FPU-trajectory's spectrum deviates from the g-breather's one at modes corresponding to a higher order 
in the excitation sequence. 

3.4.1 The original FPU trajectory: how many frequencies? 

We finally discuss a (/-breather constructed by PL series following an initial excitation resembling the one 
of the FPU-trajectory considered in the classical experiment of Fermi Pasta and Ulam in [I], that lead to 
the observation of the FPU-recurrences. 
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Figure 7: FPU-a system with a = 0.33, N = 32, E = 0.00465079 with seed mode q = 25. In panel (a) 
it is shown the evolution of the normalized instantaneous spectra E q (t) /E of modes 25, 14, 11, 28, 3, 22, 17 
for the PL solution of the (/-breather (blue) and the numerical integration (PLn) of the initial condition 
Q q L (0), P^ L (0) (orange). In panel (b) is the exponential profile of the normalized averaged energy spectra 
E q /E versus q of PL with blue spheres, of PLn with orange triangles and of the FPU-trajectory with green 
triangles. 



For the construction of the g-breather, we take a = 0.33, N = 32 and E = 0.014915. The seed mode 
is go = 1 an d truncation order of the PL series ko = 242. The chosen frequency value u)\ and amplitude 
Ai are given in HI 

Computing the normalized energies E q (t)/E, q = l,...,N — 1 for the g-breather and the FPU- 
trajectory (FiglH^a) and (b) respectively), we observe the striking difference in their time evolution. Namely, 
in the case of the FPU-trajectories, the energies E q (t) clearly exhibit a recursive behavior, leading to a 
nearly complete return to their initial values at the time t = 12000. On the contrary, the energies E q (t) 
remain practically constant for the g-breather in the same time interval. This remark is important, because 
it shows the difference in dynamics between two trajectories with quite similar normalized averaged spectra 
(FigKc)). 

This example clearly shows that the proximity of the averaged energy spectra cannot serve alone as a 
criterion for the proximity of g-breathers and FPU-trajectories in phase space. In particular, the FPU- 
trajectory of Fig J8] seems to evolve with more than one incommensurable frequencies, a fact implying that it 
lies close to a torus of dimension higher than one. By implementing the GALI method, (FigJHJd)) we have 
an indication that the numerical orbit is governed by at least five incommensurable frequencies, since G5 is 
the last of the GALI indicators found to asymptotically converge to a constant value (neglecting numerical 
fluctuations). We note, in this latter respect, that according to |14j . an excitation as in FigJH] should lead 
to the formation of a 'natural packet' of modes exhibiting a sort of internal equipartition. However, the 
prediction for the packet width yields 

a l/2 £ l/4 „ 2 _ 7 

~ 3, which is smaller than the dimension of the torus 
as suggested by the GALI indicators. This fact suggests that the dimension of the packet, as found by an 
inspection of the energy spectra and the dimension of the underlying g-torus, need not coincide. 

4 Conclusions 

In the present paper, we study the dynamical features of low-dimensional invariant objects of the FPU 
phase space called q-tori. Our main findings can be summarized as follows: 

1) We use the method of Poincare - Lindstedt (PL) series in order to compute quasi-periodic series 
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Figure 8: FPU-a with N = 32, a = 0.33, E = 0.014915 and seed mode q = 1. Panels (a) and (b) are for 
the evolution of the normalized instantaneous spectra E q (t)/E for the g-breather and the FPU-trajectory, 
respectively. In panel (c) is the exponential profile of the normalized averaged energy spectra e K of PL 
with blue spheres, the numerical integration (PLn) of the initial condition Q^ L (0), P^ L (0) with orange 
triangles, of FPU-trajectory with green triangles and of NR with black spheres. In (d) is the evolution of 
the GALI indices G 2 , G 3 , . . . , G 7 . 
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representations of trajectories approximating the motion on g-tori, up to a high order in a small parameter. 
We give some details on the way by which appropriate values of the torus frequencies are chosen in order for 
the PL method to proceed. We also discuss some numerical tests on the level of the final approximation, 
based on the use of the GALI indicator. An analysis of the convergence properties of the PL series is 
deferred to an accompanying paper 

2) Using properties of the PL series construction, we study the phenomenon of energy localization on 
g-tori. We present a theory of propagation of initial 'excitations' in the series terms. Via this theory we 
predict theoretically the form and shape of energy localization profiles for g-tori. Finally, we compare 
the latter with the energy localization profiles of orbits lying in the neighborhood of g-tori (called 'FPU- 
trajectories'). 

3) As a case of particular interest, we study low-frequency packet excitations. In this case, we show 
that g-tori solutions predict the appearance of exponential localization profiles, with a slope depending 
logarithmically on the system's parameter the a model and on fie in the /3. Furthermore, the 
g-tori exhibit invariant profiles with respect to increasing the energy E proportionally to N. Thus, the 
normalized energy spectra are invariant with respect to the re-scaled wavenumber k = q/N, as long as the 
fraction of modes excited s/N and the specific energy are kept constant. Via leading order estimates in the 
PL series, we derive a law for the slope of the exponential energy profile, both in the FPU-a and FPU-/? 
models. This turns out to be similar to that of g-breather solutions in which the 'seed mode' go is allowed 
to vary proportionally to N. In fact, such laws suggest the invariance of the g-tori's localization profiles as 
we approach towards the thermodynamic limit, provided that s/N and specific energy e are kept constant. 
However, the existence of g-tori as N — > oo is still an open issue, as we have no proof of the convergence 
of the PL series in a finite domain of initial conditions as we approach this limit. 

4) Regarding the FPU-trajectories lying in the neighborhood of g-tori solutions, we find that they 
exhibit energy localization phenomena leading to an averaged energy spectrum tending to saturate to a 
form quite close to that of a g-torus solution with similar initial excitation. We find that the process of 
energy transfer in FPU-trajectories corresponds to the so-called 'first stage' of the FPU dynamics in |22j . 
Furthermore, we provide numerical evidence that the groups of modes involved in this process are the same 
as the groups defined by the sets T>k for g-tori. A theoretical interpretation of this phenomenon remains 
an open question. 

As a final remark, we note that the time for which 'metastable states' of FPU-trajectories started close 
to g-tori persist, as well as its dependence on the system's parameters, is an interesting open question that 
can be considered as complementary to the study in [19] (where the initially excited packet is smaller than 
the natural one). This is proposed as a subject for future study. 
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APPENDIX 



A. Algorithm of determination of frequency values for the PL series construction 

We give an iterative algorithm of determination of numerical values of the frequencies u qi , qi S T>q, for 
a g-torus solution with initial excitation T>q, such that the solution constructed at every step corresponds 
to a higher specific energy than the solution constructed in the previous step. The algorithm consists of 
the following steps: 



Initialization, i) We choose some value of 'trial amplitudes' A l " al and define 

^tH A l^ ,i = l,...s. 



^Ltrial 



N ^ 2 

i=l 



(28) 



ii) We compute 'trial' frequencies by the lowest-order frequency correction terms in the series (jlOp 
corresponding to the trial amplitudes A q . tria i. 

In the FPU-a, where the frequencies oj Qi appear in the denominators of the lowest order terms, we 
implement a two-step substitution-iteration procedure. Namely we set 
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and set w^tHoi = ^ g + y?LofJ rial . 



In the FPU—/?, we simply have 



^q, trial — + ^ C^-gi ,toa£ j • • • ; Aq a ^ r i a {] 



(29) 



iii) We compute the PL series for the excitation T>q, using as frequency values Lo qi = LO qh t r ial- We 
attempt to determine numerically a root of Eqs. (|10p for the amplitudes A qi , with a root-finding technique 
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starting by A qi j r iai as guess amplitudes. If this fails, we return to substep (i), trying some lower value for 
the amplitudes A qutria i, until a successful solution is found. We store the pairs (to qi ,A qi ) of the latter, 
iv) we repeat the process (i) to (iii) for some neighboring trial amplitudes A 1 i T ^ a \ and store the pairs 



(uj' ,A' ) for the corresponding solution. We define Aq = e 



Ltrial e l,trial- 



Iteration, i) We define 



1/2 



(30) 



,i=l 



and denote by (uj qij Q, A?,,o) the original solution (u) qi ,A qi ). 

ii) We compute s neighboring solutions corresponding to the set of frequencies 0J qit j = to qi + 5ijAuj. 

iii) We compute the matrix J u = J^ 1 , where J a is a matrix defined by the finite differences 



A 2 - A 2 
Alo 



iv) Finally, we compute the next set of frequencies to be used in PL series construction by 

2NAet ^ J u 



u q% = u qi + 



(31) 



(32) 



This completes one full step of the iterative algorithm of determination of frequencies. We note that 
a change of the frequencies as in Eq, (|32p leads to an increment of the total energy corresponding to each 
successive step by an amount of order Ae^. 



in the case of the FPU-a model (see |10| for the 



B. Proof of the proposition of subsection 12.31 

We give the proof of the proposition of subsection | 
proof in the case of the FPU-/3 model). 

The proof follows by induction. Let be the set defined in Eq. (|14p . which corresponds to the set of 
all modes for which the r.h.s. of Eq. fllip is non-zero at the k— th order. According to Definition 1, the set 
of modes excited at the k-th. order is given by = \ Uo<j<fc-i w here r(k) = k + 1. For k = 1, the 
r.h.s. of Eq, (fTT|) is non-zero if B qqiq2 ^ 0, implying that for all modes q £ we have that q is either of 
the form q =| a^q^ \, or of the form q = |2iV— | (jWqW \ |, where the allowable combinations of values 
of g (2) = (qi, q 2 ) £ V\ and of ct (2) = (<Ji, cr 2 ) G S 2 , are those leading to 1 < q < N - 1. 

Assuming, now, the proposition to be true at order k — 1, one finds that, at the k-th. order, the r.h.s. 



of Eq. (fTT]) is non-zero if B q i m ^ 0, Qj; / and Qm 2> ^ 0, where n± + n 2 = k — 1. For the modes I G M, 



(na) 



and m G M n2 we then have: 
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where v; 



^ a (n i + l) q (n l + l)^ +N _ 1 

2N 



q =| I ± m | or q = \2N- 
in a combined form as: 



However, the condition B q i m ^ implies that q is necessarily of the form 
/ ± m | I . Provided that 1 < q < N — 1, the two latter equations can be written 



\2Ng- \l±m 



(33) 
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where g = or 1. Then, Eq. (|33p takes the form 

q = \2Ng± 2v x N ± 2u 2 N^ \ a {ril+l) q {ni+l) | T | (7 (n2+1) g (n2+1) | | 

= \2N(\g±u 1 ±u 2 \)-\^ k+1 k (k+1) \\, (34) 

* v ' 

!/ 

after a possible sign reversal within | • | (not affecting the absolute value) and with v a positive integer. 

However, by the restriction 1 < q < N — 1, one necessarily has that v = [(| a^ k+l ^q^ k+l ^ \ +N — l)/2iV~] . 
This concludes the proof of the proposition. 

C. The mapping K,f ) for FPU a 

We give the expressions of the quantities appearing in the mapping (|19p . up to the order k = 3. We 
have 
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where go, gi, ^2 and q^ represent any mode belonging in the sets Vq, T>%, T>2 and V3 respectively. 
D. Explicit expressions for the leading order terms of the PL series 

(k) 

We prove by induction that the leading order terms Q q for the modes q £ 2\- are given by Eqs. (H 
with the quantities TZ q h \ )C q k ^ given by (fl~8|) . (fl~9j) and p0|) . We focus again on the FPU-a model. 
For k = 1, the solutions of Eqs. (jlip read 
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(35) 

so, Q q l) satisfies Eq.(TT7I). 

Assume now that Eq. (|17|) holds true for the solution Q q (t) at the order k — 1. For simplicity, we use 

the notation n\ ni) = TZ\ nz) (n^ n ^) and /c[ ni) = /cj ni) (n( r(ni)) ), i = 1, . . . ,r(fc). At order fe, Eq. (HH) takes 
the form 

Q<£)+nlQ(f)= (36) 
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By replacing the term 
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into the above equation, one has that 
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However, the solution of the latter equation is 



<4 fe) (*) 



Thus, Eq. (fT7|) holds true at the order k, with the expressions 1Z q k \ K, q given by (fT8|) . (fT9|) and ([20 
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E. Localization profiles of q tori estimated by leading order terms in the PL series 

We derive estimates for the form of the energy localization profiles for g-tori solutions corresponding 
to an initial excitation of the modes 1 < q < s, with s varying proportionally to N. 
We first make the following estimates: 



i) For any mode q £ T>^ we use the approximation q ~ c^s, where c^s is the mid mode of i.e. 
q = ks+ [s/2] (c fc ~ A; + 1/2) in FPU-a and q = 2ks (c k = 2k) in FPU-/3. 

ii) For the unperturbed frequences we use the approximation £l q ~ irq/N. 

hi) For m £ T>k and for 'almost resonant terms', for which m — av( k '>n = 0, we use the approximation 
|n m - <J (rik)) uJn {k)) | ^ vr 3 m 3 /(24iV 3 ) (See Appendix B of [ID] for its derivation). 

iv) For the quantities &m of Eq. ([2T]) we set 

|£(*)(n(K*)h| ~ 9m 12iy2 

v) Finally, for the total energy of the system we use the approximation E ~ Y2n=i s^ n i wnere 
E n ~ 1/2^0*. 

We will now derive an estimate for the quantity \Q q of Eq. (|17p . We first write some approximations 
for the terms of Eqs. (|18|) . (|19|) and (|20|) . Recalling that r(fe) = fc + 1 in the a model, and r(k) = 2k + 1 in 
(3, and taking into account the approximations (i)-(v), we find 
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The quantities fZ q k a and /C^ of (fl~9j) and (f2"0j) cannot be evaluated analytically. However, as already 

mentioned in section 12.41 their form is a polynomial of degree k — 1 in the terms £m , I = 0, 1, . . . and a 
polynomial of degree k in B q ^ m (or Cqj^n). We denote here by "P(<7, ra( r \ m( fc_1 )) the product of k factors 

of the coefficients B q ^ m (or C q)l , m>n ) in K$l (or /C^). 
The size of the mid mode is then 
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where by we note constants in powers of A;. Due to the term "P((7, n^'\ m^ -1 )), the sums over 
and m^ _1 ) in (14=11) give rise to a factor s^ r ^ +1 ^ 2 , i.e. s fc / 2+1 in a and to in j3 @. Replacing mj, 
i = 1, . . . , k — 1 and g in (|4ip by their mid mode expression c^s, one has 

// N 2k+1 r(fc )/2 / A „s 



The energy in each group of modes is then estimated as 

1 94-^9.9 1 Oh /vrc fc s\2 /„„ j\T\ 2 ( 2fc + 1 ) 

7TS / 

72 \ 2fe 



r«0 



C: • ( ^1 ) . (43) 



Vr 2 ^ 2 



Finally replacing r(k) and \i by and aj\JlN for FPU-a, or 2/c+l and /3/(2N) for FPU-/3 respectively, 
we arrive at the exponential laws 

a 2 EN 3 \ k /a 2 eN 4 \ k 



F. Precise values of frequencies and amplitudes in all paper's numerical examples 

4 These factors are found by recursively solving ^2 n ( r ) GT> r V{q, n^ r \ m' fe_1) ), for those q and m'*" 1 ' that maximize the 
results. In /3 one finds for k = 1 the factor 2s(s — 1), for k — 2 the factor 4s 2 (s — 1), etc. 
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Example 



Frequencies 



Amplitudes 



Fig.l, 2? = {1,2,3,4} 

Fig.4a, £> = {1,11, 21, 31} 

Fig.4b, V Q = {60, 61, 62, 63} 
Fig.4c, X>o = {63,64,65} 

Fig.4d, V Q = {94, . . . , 98} 

Fig.6, A, = {1} 
Fig.7, V = {25} 
Fig.8, V = {1} 



wi = 0.09813690483108113 
w 2 = 0.19603499627351076 
lj 3 = 0.29346131418681830 
lo 4 = 0.39018079936090050 

wi = 0.09814720863898445 
wn = 1.0282023436963477 
w 2 i = 1.7154505429645952 
w 3 i = 1.9975829222954775 

Lo 6Q = 1.9903455323483594 
w 6 i = 1.9945569202421924 
w 62 = 1.9975668544869103 
w 63 = 1.9993735255838023 

w 63 = 1.3967520180075357 
w 64 = 1.4142130013133483 
w 65 = 1.4314611496716776 

w 94 = 1.8284112766464593 
w 95 = 1.8382194182395468 
w 96 = 1.8477507307792000 
u g7 = 1.8570037790005570 
w 98 = 1.8659771693147402 

ui = 0.09814109959448596 

w 25 = 1.8830741847088537 

oj! = 0.09844127049688513 



A 1 = 0.098097447281306230 
A 2 = 0.048959912346992665 
A 3 = 0.032444239210574590 
A 4 = 0.023892612612360370 

At = 0.2699420271816366000 
An = 0.027437337334704704 
A 21 = 0.015901830855243902 
A 31 = 0.014550407721151376 



-460 
-461 
A 6 2 

-4 63 

-4 63 
A 6A 
A 65 

A m 
A 95 
A g6 

A97 
A qR 



0.010626655362545707 
0.010540492500677328 
0.010486090307194270 
0.010462240436802113 

0.005128073559888574 
0.005045052412361571 
0.004963639646819611 

0.008845316872292220 
0.008752766020034742 
0.008700391793859069 
0.008662026577864702 
0.008661383820872627 



Ai = 0.18574105489382606 
A 25 = 0.0512142216322969 
A x = 1.3009149083500795 
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